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Bright sink-type localized states in exciton-polariton condensates 
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The family of one-dimensional localized solutions to dissipative nonlinear equations includes a 
variety of objects such as sources, sinks, shocks (kinks), and pulses. These states are in general 
accompanied by nontrivial density currents, which are not necessarily related to the movement of 
the object itself. We investigate the existence and physical properties of sink-type solutions in non- 
resonantly pumped exciton-polariton condensates modeled by an open-dissipative Gross-Pitaevskii 
equation. While sinks possess density profiles similar to bright solitons, they are qualitatively differ¬ 
ent objects as they exist in the case of repulsive interactions and represent a heteroclinic solution. 
We show that sinks can be created in realistic systems with appropriately designed pumping pro¬ 
files. We also consider the possibility of creating sinks in a two-dimensional configuration with a 
ring-shaped pumping profile. 

PACS numbers: 71.36.-l-c, 03.75.Lm, 42.65.Tg, 78.67.-n 


I. INTRODUCTION 

Strong coupling of semiconductor excitons to micro- 
cavity photons results in the appearance of spectral res¬ 
onances associated with mixed quantum quasiparticles 
called exciton-polaritons^. These particles exhibit ex¬ 
tremely light effective masses, few orders of magnitude 
smaller than the mass of electron, which allows for the ob¬ 
servation of physical phenomena related to Bose-Einstein 
condensation already at room temperatures^”®. At the 
same time, polaritons exhibit strong exciton-mediated 
interparticle interactions and picosecond lifetime due to 
their photonic component. They are actively studied 
both from the point of view of fundamental interest®”® 
and potential applications 

In recent years, great attention has been devoted to the 
study of nonlinear self-localized states of superfluid po¬ 
laritons, such as dark and bright solitonsSolitons 
are nonlinear wavepackets which preserve their shape 
thanks to the balance between dispersion and nonlin¬ 
earity^®. They have been applied to long-distance op¬ 
tical fiber communication^® as well as description of nu¬ 
merous physical systems. Polariton solitons have been 
demonstrated both in the cases of resonant^®d3 
resonant pumping^'^”^^. To date, no bright states were 
shown to exist in the nonresonant case with homogeneous 
pumping. 

Polariton superffuids are inherently nonequilibrium 
systems in which the balance between pumping and loss 
is an essential factor®’®’®®. In many of the previous stud¬ 
ies, this aspect was treated as an unwanted complication 
of the theory. Standard models, such as the conserva¬ 
tive Gross-Pitaevskii equation, were frequently used to 
describe solitons. However, it is well known that self- 
localized solutions in dissipative systems have qualita¬ 
tively different properties than their conservative coun¬ 
terparts. In the case of repulsive interactions, only one 
type of one-dimensional solution exists in the conserva¬ 
tive theory - dark or bright solitons, depending on the 
sign of the effective mass. In the dissipative case, a fam¬ 
ily of qualitatively different localized states exists, in¬ 


cluding sources, sinks, shocks and pulses®^”®®. In gen¬ 
eral, they exhibit nontrivial internal density currents and 
may undergo complicated, sometimes even chaotic dy¬ 
namics'^®’®®. 

In this paper, we demonstrate the existence and 
stability of bright self-localized solutions of the open- 
dissipative polariton model. These solutions are classi¬ 
fied as sinks (anti-dark solitons)®®’®®, the name reflecting 
that their structure corresponds to terminating lines of 
incoming density currents, with a local increase of loss. 
While sink-type solutions possess density profiles similar 
to bright solitons, they are qualitatively different objects. 
In contrast to bright solitons, they exist in the case of 
repulsive interactions and represent a heteroclinic solu¬ 
tion connecting two counterpropagating plane waves. We 
demonstrate the dynamics of sink formation and their 
stability in a realistic model with appropriately chosen 
pumping profile. We investigate systematically the prop¬ 
erties of sinks and provide an approximate analytical for¬ 
mula for their shape. In the two-dimensional case, we 
show that sink creation is hindered by the spontaneous 
proliferation of vortices, which destroy the supercurrents 
necessary for the existence of a symmetric sink solution. 


II. MODEL 


We consider a polariton condensate in the one¬ 
dimensional (ID) setting, e.g. trapped in a microwire®^. 
We model the system with the generalized open- 
dissipative Gross-Pitaevskii equation for the condensate 
wavefunction ip{x, t) coupled to the rate equation for the 
polariton reservoir density, ni?(a;, t)'®’®®’®® 
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where P{x) is the exciton creation rate determined by the 
pumping profile, m* is the effective mass of lower polari- 
tons, 7c and 7 r are the polariton and exciton loss rates, 
and ^)/ V 27 rd 2 are the rates of 

stimulated scattering into the condensate and the inter¬ 
action coefficients, rescaled in the one-dimensional case. 
Here, we assumed a Gaussian transverse profile of IV'P 
and riR of width d. In the case of a one-dimensional 
microwire^', the profile width d is of the order of the 
microwire thickness. We also introduced D = 1 — iA 
with A being a small constant accounting for the energy 
relaxation in the condensate^^^^*^. 

To obtain a system of dimensionless evolution equa¬ 
tions, it is possible to rescale time, space, wavefunction 
amplitude and material coefficients according to t = rt, 
X = ^x, iIj = nR = {^/3)-^nR, = 

{^I3/t)R, {g}P,g]P) = {f>^P/T){gc,gR), [i c^ir) = 
'r~'^{lC,lR),P{x) = (l/C/3T)P(a;), where^ = y/fvrj2^, 
while T and (3 are arbitrary scaling parameters. We 
rewrite the above equation in the dimensionless form (we 
omit the tildes for convenience) 
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In the above transformations the norms of both fields 
NjIj = J li^l'^dx and Nr = J URdx are multiplied by the 
factor of j3. 


III. SINK-TYPE SOLUTIONS 
A. Description of sinks solutions 


The structure of sink solutions can be understood most 
easily as a result of interaction of two counterpropagating 
nonlinear waves. In the linear case, two waves emitted 
by distant sources give rise to a standard interference 
pattern, as shown in Fig. 1(a). In the case of a dissi¬ 
pative model with nonlinear gain or loss coefficients, the 
interference pattern can be replaced by a localized den¬ 
sity peak, sometimes exhibiting oscillating features, as 
depicted in Fig. 1(b). Sink is a heteroclinic solution con¬ 
necting two plane waves emitted by the sources at each 
side. The two waves collide at the sink position, where 
the incoming density currents are dissipated^^. 

The sink density pattern is a result of the ability of 
the dissipative medium to smoothen out density “dips” 
and “peaks”, which are present in the standard interfer¬ 
ence pattern. When one of the incoming waves reaches 
the area occupied by the other, the resulting interference 
leads to decay of waves. Let us consider one of the sim¬ 
plest dissipative nonlinear wave equations, the complex 




FIG. 1: Density patterns created with counterpropagating 
waves, (a) Interference pattern in the linear regime, (b) Sta¬ 
tionary sink-type localized pattern in a model with nonlinear 
dissipation. 


Ginzburg-Landau equation^^’^^ (GGLE) 
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This equation can be obtained from the system of equa¬ 
tions (2) in the limit of fast relaxation time of the reser¬ 
voir, which is “slaved” by the slower ^ dynamics, and in 
the linearized approximation Ri ng = ImH/ImC = 
(P — Pth)He Ho is the dynamical equilibrium density 
when the loss and gain are balanced). Here, we as¬ 
sumed a homogeneous pumping P{x) = const > Pth 
and introduced the threshold power Pth = 'Yr'IC / R- The 
GGLE parameters are B = ijcH ~ Hr + *P/2)(1 -I- 
Rno/'yA)P/lA and C = {gR + iRl2)PRH\ “ 9C with 
7 A = Png + 7fl- It is clear that the existence and stabil¬ 
ity of the homogeneous steady state with ng > 0 (which 
in general can be also a plane wave solution) requires 
ImP > 0 and ImC > 0. Under these conditions, per¬ 
turbations of the steady state with density ng exponen¬ 
tially decay, which is the reason for the above mentioned 
smoothening. Any areas of density higher than ng corre¬ 
spond to net loss, and those with density lower than ng to 
net gain in Eq. (3). Nevertheless, nontrivial (non-plane 
wave) stationary solutions can still exist^®’^^’^^. 

One can formulate another necessary condition for sta¬ 
bility of sinks based on the modulational (or Benjamin- 
Eeir) stability of the incoming plane waves^^. In the case 
of GGLE with real P > 0, this is assured by the condi¬ 
tion Re C < 0, which in terms of Eq. (2) translates into 
P/Pth > Hc9rHr9c)^ as shown recently in^^. We note 
that this is a necessary condition for stability, and an ac¬ 
tual domain of stability of plane waves may be smaller^®. 

It may seem natural to treat sinks as dissipative ana¬ 
logues of bright solitons. These states are, however, 
qualitatively different from each other. Bright solitons 
exist in the conservative limit of the GGLE (3) with 
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FIG. 2: Patterns created by polariton waves emitted by high- 
intensity sources on the two sides with pumping profiles as in 
(4) (dashed lines). Panel (a) presents a stationary interference 
pattern obtained by integration of ( 2 ) in time for parameters 
A = 0.1, R = 0.96, gc ~ 0.63, qr = 1.91, 7 c = 0.9, 7 ^ = 
0.6 at low pumping powers, (b) With a stronger pumping 
Po a nonlinear sink-type solution is created. Corresponding 
parameters in physical units are: time unit r = 7 ^^ = 3ps, 
length unit ^ = 1.9 /rm, g — 3.9 geVgrr?, R — 9 x 10“® gm^ 
ps“^ for d = 2 gm, m* = 5 x 10~®me, and /3 = 0.003. 


Im ((7,5,11) = 0, which is the celebrated Nonlinear 
Schrodinger equation^' (or Gross-Pitaevskii equation in 
the context of degenerate bosons'^^). Sinks require non¬ 
zero incoming currents for their existence^^, and exist 
in the case of a stable background. Re (CD) < 0, while 
bright solitons exist only in the self-focusing (modula- 
tionally unstable) case with CD > 0. 


B. Sinks in the exciton-polariton model 

To create sink solutions in the model described by (2), 
one has to provide sources of counterpropagating waves 
as described above. We consider the following pumping 
prohle created by a pumping beam with spatially varying 


intensity 

P(X) = - (Pmax " (4) 

as shown in Fig. 2 (dashed lines), where we used the 
smoothness parameters a = 100 and /? = 80. The profile 
exhibits jumps of P(x) at x = iLwb and x = ±Ws, typical 
of the Super-Gaussian terms in (4). 

The sink is created in the central area with pump¬ 
ing intensity P = Pq. The side areas with P = Pmax 
are the sources of polariton waves. This flow is obtained 
thanks to the repulsive polariton-polariton and reservoir- 
polariton interactions QCt Qr > 0, which create an effec¬ 
tive potential hill in the areas with high pumping density 
P = Pmax- In result of interaction of the two nonlinear 
waves propagating towards the center, depending on the 
parameters of the system, interference or localized sink 
patterns can appear, as shown in Fig. 2(a) and 2(b), re¬ 
spectively. Note that at the interface between the high- 
and low-pumping areas at a; = ztwg in 2(b), oscillating 
time-dependent states form, which do not, however, pre¬ 
clude the formation of sinks. We obtained also different, 
more complicated nonlinear patterns for other values of 
parameters, especially in the case when the two sources 
were relatively close to each other. These patterns did 
not have a localized character such as the one shown in 
Fig. 2(b). A small value of relaxation coefficient A = 0.1 
was in some cases necessary to attenuate high-momentum 
modes in simulations and obtain physically relevant so¬ 
lutions. 

Figure 3(a) shows the dynamics of the sink creation 
process. Initially, we assumed zero density of excitons 
and a small white noise in the polariton field, but we 
checked that the final stationary state is practically inde¬ 
pendent of the form of the initial condition. First, in the 
area between the sources a condensate is created with ap¬ 
proximately zero momentum. The wave fronts generated 
by the sources gradually move towards the center, where 
they collide creating the stable sink structure. The waves 
are “stopped” by the sink due to the nonlinear character 
of the gain and dissipation. 

It is important to note that the sinks are in general 
not completely stationary, as in Fig. 3(a), but may be 
put in motion by the imbalance of momenta of the waves 
emitted by the two sources, see Fig. 3(b). In this case the 
sink moves with a constant velocity proportional to the 
mismatch between the two wavevectors^^. If the balance 
is restored after some period of time, the sink stops at 
the new position. The sink may be also stopped after 
reaching the weaker source, as demonstrated in Fig. 3(b). 

IV. DOMAIN OF EXISTENCE AND 
PROPERTIES OF SINKS 

In this section we describe a systematic investigation of 
stationary sink solutions of the exciton-polariton model 
with homogeneous pumping, P(x) = const. We sub¬ 
stitute tj}(x,t) = 4’{x)e~''^* and nji(x,t) = nR(x) into 
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FIG. 3: (a) Evolution towards a stationary sink state from a 
small initial noise. Parameters are as in Fig. 2(b). (b) The 
case of asymmetric pumping profile, with momentum mis¬ 
match between the waves from the two sources. The sink 
moves towards the weaker source, where it is stopped but not 
destroyed. 


Eqs. (2) to obtain a single ordinary differential equation 
for the profile of the stationary state 
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We complement the above equation with boundary 
conditions. At x = 0, which we chose to be the sym¬ 
metry point without loss of generality, the first deriva¬ 
tive is equal to zero, d())/dx = 0. At x = -boo the so¬ 
lution tends to the plane wave with the norm equal to 
|(()p = no = {P — Pth)he (hr practice we impose this 
condition on the last point of the computational mesh). 

We solve the boundary problem with the shooting 
method using the Newton minimization algorithm. We 
keep ^|a;=o = 0 and change the value of ^(0) while solv¬ 
ing (5) with the Runge-Kutta algorithm on a certain in¬ 
terval 0 < X < Xniax- The Newton method is then used to 
find a solution that satisfies boundary condition = ng 
at X = Xmax within a given tolerance. We then extend 
the interval boundary x^ax slightly and repeat our pro¬ 
cedure. This method proved to be an effective way to 
obtain a localized state on a large x domain. 



FIG. 4: Phase diagram depicting parameters for which a sink 
solntion was obtained with (5) (gray area). Parameters are 
R = 0.76, gc = 0.38, qr = 0.76, 7 c = 0.75, 7 n = 1. 



FIG. 5: Panels show the values of |'0min max (x)p, XI 

and the wave vector k (see Fig. 1 for descriptions). On panel 
(a) we keep a constant value of P = 1.45 but vary the chemi¬ 
cal potential g. The value of Xmin quickly decreases indicating 
that at higher the sink is more densely undulating. The de¬ 
pendence of |'i/’min(x)p and |'!/)max(x)|^ shows that the sink 
“height” is increased at higher g. In panel (b) we change the 
value of P while /r is chosen slightly higher than the lower 
threshold fj-min for which a sink solution occurs. The depen¬ 
dence of |'i/)inin(x)|^, |V’max(x)|^ On P indicates that sink height 
is approximately constant while the background density no 
grows with P. 
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Figure 4 presents a phase diagram showing the domain 
of existence of sink-type solutions in the parameter space 
of P and /i with values of other parameters fixed. The 
shaded area, corresponding to parameters for which the 
algorithm converged to a sink solution, is limited from 
below by the natural minimum given by the chemical 
potential of the steady state 

_ 2 , PgR ia\ 

Mmin — gCR-O + / i p 2 \ ’ 

{jR + Rn^) 

The range of /r for which the sink solution exists turned 
out to be the largest for moderate pumping intensities 
P Ri 1.5Pth- Although the range of /r appears to be small, 
it corresponds in fact to a broad range of wavenumbers 
of incoming waves, from approximately zero to around 
6, as shown in Fig. 5(a) with a dash-dotted line. In 
this figure, the vertical line corresponds to the minimal 
value of fi, given by Eq. (6). The dependence between 
the chemical potential and the wavenumber of incom¬ 
ing waves can be calculated by taking the x —>■ ±oo 
limit away from the sink, for which Eq. (5) reduces to 
^ = k'^ + gcRo + PgrUir + R'rV) under the condition of 
balanced gain and loss. It is then clear that the neces¬ 
sary condition for the existence of sinks is that the kinetic 
energy is much smaller than the nonlinear energy of the 
wave. 

The dashed and dotted lines in Figure 5(a) depict the 
minimum and maximum density of the sink profile, see 
Fig. I. With increasing /i, which corresponds to increas¬ 
ing k vectors of incoming waves, the sinks become larger 
and more highly modulated, while in the g, = ^min limit 
they transform smoothly to the flat homogeneous state 
4>{x,t) = The position of the first minimum 

of the density Xmin (solid line) decreases with the increase 
of fjL, which is related to the increasingly dense interfer¬ 
ence pattern of the tails of the sink. 

Figure 5(b) shows the dependence of the same proper¬ 
ties of the sink as described above, but with increasing 
pumping power P. Here g, is kept at a slightly increased 
level with respect to g-min- This corresponds to a con¬ 
stant small value of k. In this case, the increase of P 
leads to the increase of average sink density, but with¬ 
out almost any change of the difference l^imaxP — 

On the other hand, the value of Xmin shows a strongly 
nonmonotonous character, decreasing for small P and in¬ 
creasing again for large P. This shows that in the case 
of small k the position of first minimum is not simply 
related to the incoming wavevector k. 


V. ANALYTICAL SOLUTION 

If the “height” of the sink is relatively small compared 
to the steady state density, | Itpix, t)\—no \ ng, and the 
amplitude {ipix^ t)| is slowly varying in space, an approx¬ 
imate analytical solution can be found^^. We rewrite the 
steady-state solution in a homogeneously pumped con- 
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FIG. 6: Comparison of the numerical sink solution from the 
shooting method with the approximate analytical solution for 
two sets of parameters. Parameters are R = 0.76, g = 0.38, 
qr = 2g, 7c = 0.75, 7 _r = 1, P = 2.0, fi = 1.2683 for (a) and 
P = 1.85, g = 1.2173 for (b). 


densate applying the Madelung transformation for '0(x, t) 
nR{x,t) = URix), 


where a{x) is the amplitude, ip(x) is the phase and g is 
the chemical potential of the condensate. Neglecting the 
spatial derivatives of a{x), Eq. (2) can be rewritten as 


ig 


ga^{x) P gRnR{x) -b 
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- {RnR{x) - 7c) 



( 8 ) 


The real part of Eq.(8), which can be interpreted as the 
continuity equation, gives 


nR{x) 


2 

R 



ic 
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On the other hand, solving the equation for reservoir den¬ 
sity (2) in the steady state we get 


nR{x) 


P 

a?{x)R + ^R 
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We expand this formula into Taylor series of degree two 
around a^i^x) = ug = {P/lc) - {ir/R) 
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Comparing Eqs. (9) and (11) we obtain 


, - PRlC + ICIR 

a (X) =-- (12) 

From the imaginary part of Eq. (8), using (9) and (12) 
we get 





PRlc + iclR 
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With the definition ^(x) = ^(/^(x) we obtain the first 
order differential equation 
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With the solution 


i?7c \ 2 a , 


(16) 


where S = R{PRgjc - Rgic - giclR + gRlc) a = 
gPR — gRJc- Calculating the amplitude with (12) we 
finally obtain 


a^(x) 
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1 + tan^ 


PRjC + 7c7ii] 


/ 1 7cV^ 
1^2 a 


+ 

(17) 


Comparison between the analytical and numerical so¬ 
lution obtained with the shooting method is shown in 
Fig. 6. In general, very good agreement is obtained for 
small fc, when the sink profile is flat and broad, with 
no oscillating tails. The tails obviously cannot be repro¬ 
duced by the approximate solution (17), which is visible 
in the (b) panel of Fig. 6. 
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FIG. 7: Snapshot of a two-dimensional solution of (19) after 
a long time of evolution t = 1600, generated by a pumping 
profile in the shape of a ring. Quantum vortices are clearly 
visible in the density (top) and phase (bottom) of 
Sink solutions in 2D are absent due to proliferation of vortices. 
Parameters are A = 0.1, R = 0.96, gc = 0.63, qr = 1.91, 
7c = 1, 7fl = 0.6 with Po = 6, Pmax = 10,ri = 155, r 2 = 220. 


VI. TWO-DIMENSIONAL CASE 


Additionally, we performed a series of simulations to 
investigate whether sink creation is possible in the two- 
dimensional (2D) version of the exciton-polariton model, 
described by the equations 



/ \ 
yclx^ dy'^j 


+ gc\'>P\’^ + gRnR+ 


+ 2 “ 9'c) 


4’, 


duR 

dt 


= P{x, y) - {jR + R\4j\'^) riR, 


(18) 

(19) 


where gR,c, ^R,C, R, and P{x,y) are dimensionless pa¬ 
rameters obtained from physical ones in an analogous 
way as in the ID case. 

One of possible choices of the pumping profile corre¬ 
sponds to a constant P = Pq pumping intensity on a 
circle of radius ri, surrounded by a ring of P = Pmax 
with inner and outer radius xi and r 2 , respectively. In 
Figure 7 we show a typical state obtained with this kind 
of pumping profile and an evolution from a small ini¬ 
tial noise. In general, after a certain time of evolution, 
smaller or lager number of vortices are spontaneously cre¬ 
ated, and often a stationary state could never be reached 
even with a very long integration time. Vortices may 
appear spontaneously during condensation in a process 
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analogous to the Kibble-Zurek mechanism^"^"^® as well as 
due to the emergence of supercurrents in an inhomoge¬ 
neous system^’^'. Despite using various combinations of 
system parameters, as well as asymmetric ring pumping 
profiles, we were not able to obtain any stable structures 
that would resemble one-dimensional sinks of the previ¬ 
ous sections. 

We note that similar pumping profiles were used in 
several experiments^®’^® where multi-lobe or vortex pat¬ 
terns were observed. However, the experimental patterns 
were in most cases regular, which suggests that they cor¬ 
respond to the linear regime as in Fig. 2(a). In the 
case of strong nonlinear interactions, regular vortex chain 
patterns could be observed with the resonant pumping 
scheme®®. Depending on the system parameters, these 
could be destroyed by spontaneously nucleating vortices 
created through a hydrodynamic instability, which is con¬ 
sistent with our simulations. 


VII. CONCLUSIONS 

We demonstrated the existence and stability of a fam¬ 
ily of bright sink solutions of the open-dissipative polari- 


ton model. In contrast to bright solitons of conserva¬ 
tive models, sinks exist in the case of repulsive interac¬ 
tions and are created in a collision of counter-propagating 
waves. We studied the dynamics of sink formation in a 
realistic one-dimensional polariton model with appropri¬ 
ately chosen pumping profile. We studied the domain of 
existence of sinks in parameter space and their physical 
properties. An approximate analytical formula for the 
sink shape, valid in the case where sinks do not possess 
oscillating tails, was determined. In the two-dimensional 
ring-shaped confguration, sink solutions were not found 
due to the spontaneous appearance of vortices. 
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